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Characterization and quantification of multipartite entanglement is one of the challenges in state- 
of-the-art experiments in quantum information processing. According to theory, this is achieved via 
entanglement monotones, that is, functions that do not increase under stochastic local operations 
and classical communication (SLOCC). Typically such monotones include the wave function and 
its time-reversal (anti-linear operator formalism) or they are based on not completely positive maps 
(e.g. partial transpose). Therefore, they are not directly accessible to experimental observations. We 
show how entanglement monotones derived from polynomial local SL- invariants can be re- written in 
terms of expectation values of observables. Consequently, the amount of entanglement - of specific 
SLOCC classes - in a given state can be extracted from the measurement of correlation functions 
of local operators. 



Introduction. - The preparation of multipartite quan- 
tum states and the characterization of their entangle- 
ment properties is one of the central issues of ongoing 
experimental research in quantum-information process- 
ing. Technological progress has made it possible to gen- 
erate, e.g., GHZ and W states for three and more qubits 
with trapped ions (up to eight qubits, [U, 0) a-nd pho- 
tons (up to six qubits, 0, Q). Moreover, an entangled 
four-qubit state that is locally inequivalent from the GHZ 
state, the so-called cluster state could be generated with 
photons 0, Q. 

The experimental detection of entanglement is a rather 
subtle issue 0, Until now, the experimental char- 
acterization of entanglement has been carried out using 
entanglement witnesses or Bell inequalities (for a recent 
review see Q). The advantage of these methods is that 
the presence of entanglement in a given state can be in- 
ferred simply by measuring certain correlation functions 
of local observables. Although these methods are not 
strictly quantitative since to date there is no formulation 
of these concepts such that their results arc invariant un- 
der stochastic local operations and classical communica- 
tion (SLOGG) [l^, recent proposals exist on how (tight) 
lower bounds can be extracted with variable numerical 
effort ^M^- 

On the other hand, necessary criteria for the quantifi- 
cation of entanglement have been worked out and lead to 
the concept of the entanglement monotone [l^ . Essential 
properties are local SU{2) invariancc, convexity on the 
space of density matrices and non-increasing value under 
local operations and classical communication (LOGG). 
The non-increasing property needs to be extended to in- 
clude also stochastic local operations and classical com- 
munication (SLOGG) [l^, [ig. The connection between 
SLOGG and invariance under local SL{2, C) has been 
established p^l-[l7j. Based on this, the two SLOGG- 
inequivalent entanglement classes for three qubits have 
been discovered whose representatives are the W and the 



GHZ state [l5[. While the former state exclusively con- 
tains pairwise entanglement (measured by the concur- 
rence [131 ) , a measure for the latter has been derived for 
pure states, the 3-tangle[l9|. Subsequently much work 
has been devoted to the classification of iV-partitc qubit 
entanglement {N > 3) (20l - [26j . but up to now it is not 
clear what the essential criteria for such a characteriza- 
tion might be. However, it is known that local SL{2,C) 
invariance plays a major role [l^l, and 5'L(2, C) invari- 
ants - both from classical invariant theory [23, |28| - |32| as 
well as from other methods [13, - can be written in 
terms of the coefficients of a state alone (and not their 
complex conjugate in contrast to SU(2) invariants) and 
thus are naturally represented as expectation values of 
certain antilincar operators [l^, ■ This has the draw- 
back that those quantities are not directly accessible in 
an experiment. Interestingly, also other prominent en- 
tanglement measures related to not completely positive 
operators (as the Percs-Horodecki criterion [l^jllal) share 
this draw-back, which might indicate that sensible mea- 
sures for entanglement are hard to access experimentally 
for conceptual reasons Q, [s^ . For measures as the neg- 
ativity the obstacles from dealing with not completely 
positive maps can be circumvented by adding the iden- 
tity map [331 or by certain contraction techniques [38| . 
The addition of an identity map has also been applied to 
the concurrence [39|. In this work, we demonstrate how 
the antilinear formalism in Ref. |2 3| . | 2 4ll - and hence all 
polynomial S'L(2,C)®* invariants [SJl - can be directly 
rewritten in terms of linear operators and expectation 
values thereof by means of an invertible mapping j40l |. 
This combines the advantages of experimental accessi- 
bility and controlled quantification of entanglement (in- 
timately related to SL invariance). It provides experi- 
mentalists with a tool to directly measure genuine mul- 
tipartite qubit entanglement. The mapping itself does 
not require qubits as local entities and therefore it could 
readily be applied to future SL invariants also for higher 
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local dimension (e.g. for higher spin). On the other hand 
a connection to alternative approaches based directly on 
expectation values of linear operators (see e.g. [4l|) can 
be established by means of the inverse mapping. This 
comparison also highlights the advantage of an antilin- 
ear formulation as far as the construction of e.g. com- 
plete sets of invariants is concerned for future progress in 
entanglement classification. 

The paper is organized as follows. After briefly review- 
ing the antilinear formalism we present the connection 



J 



with linear expectation values. Then we discuss some 
basic conditions to measure such expectation values in 
an experiment. 

Antilinear formalism. - In Refs. (23l.[23| we have shown 
that SLOCC invariants can be obtained via the expecta- 
tion values of antilinear operators, as a straightforward 
generalization of the well-known entanglement measures 
concurrence C(0) = |(0|o-y (g) cry£|0) | for a two-qubit 
state 10) and 3-tangle r3(x) for a three-qubit state |x): 



T3(x) = \{x\<yf.(^Ty(g,ay(^\x){xW''(E>(Ty(g,ay€\x)\ 

= Kxl • (XIK ®<^y® • K ®CFy® '7y)<^lx) • lx)l 
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Here the symbol • indicates a tensor product where the 
jth factor is related to the jth copy of \x), cr^i are the 
Pauli matrices {ay = (T2, ctq = 1) and £ is the complex 
conjugation operator. The "metric" for the contraction 
of lower and upper indices is g^j^u = diag {1,-1,0,-1}. In 
Eq. ([2]) it is evident how the 3-tangle can be written as an 
expectation value of an antilinearly hermitian operator 
with respect to a multiple copy of the state |x)- The 
basic ingredient for the construction of global antilinear 
operators ("filters") are local operators ("combs") that 
have zero expectation value for any pure state of the local 
Hilbert space. The two simplest independent combs are 



Oi = ay€. 



(3) 



O2 = o-^ • cr^'e: := ^ .gp^cTp • (T^c: . (4) 



Linear expectation values.- Consider now an g-qubit 
state I ■(/')• Any polynomial g-qubit SL®'^ invariant (an 
antilinear expectation value of several copies of the pure 
state state) can be expressed in terms of linear observ- 
ables after rewriting the product of the (complex) g-qubit 
invariant and its complex conjugate expression 
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Here we have rearranged the factors and introduced the 
symbol o to indicate the tensor product of the jth copy 
and the corresponding copy from the complex conjugate 
expression. The square brackets enclose - for the kth 
qubit - as many pair-wise o products as there are copies 
of the state lip) in the antilinear expression of the invari- 
ant. The permutation operator Pj,fe exchanges the jth 
copy of the kth qubit with the jth copy of the kth qubit 
in the complex conjugate expression [49[ and the product 



extends (for the kth qubit) over all copies j of that qubit. 

Now, even the complex conjugation in the copies \ip*) 
can be removed by including a partial complex conju- 
gation acting on all operator factors to the right of a o 
symbol. It is worth mentioning that since we deal with 
Hermitean operators, this partial complex conjugation is 
equivalent to a partial transpose. As an illustration of 
the presented invertible mapping we quote the result for 
the concurrence 
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*){(t)*\(Ty ® (7y 



M„ 



diag {1, -1, -1, -1} 



(6) 



In the above equation, the symbol "o" stands for a usual 
tensor product; we use a different symbol for highlighting 
that different copies of the same state are multiplicatively 
combined where the different copies arise from taking the 
absolute square of an SL invariant. In contrast, "•" indi- 
cates tensor products of different copies of the pure state 
leading to a given homogeneous degree of the (complex) 
polynomial SL invariant (see above), and "(8)" represents 
the tensor product structure of the multi-qubit Hilbert 
space into single qubits. Note that the metric M^^ for the 
summation in Eq. ([6|) is precisely the Minkowski metric. 

Hence, we have found a rule to write the linear operator 
corresponding to a comb O ~ O € as 



£[6] := {to€){dod)I' 



(7) 



where the complex conjugation acts only on the operator 
factors in O and P to the right of the o symbol. For the 



comb (T^ • (T^ we find 



(8) 



where we have defined 



©kA^iv — '^K+A.S'^A+i/.skKM^Ai/l^^i/ + S^fi^Xi'Jtj.L' (9) 
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These results represent a prescription to express arbitrary 
polynomial SL{2, C)®'' invariants in terms of expectation 
values of linear operators. The principal idea of the map- 
ping ([7|) is generic to SL invariants and not limited to 
qubits. Clearly, the elegance of the antilinear formalism 
can be noted by comparing Eq. ([T]) to the more cum- 
bersome representation of the 3-tangle in terms of spin 
correlation functions (the formal existence of such an ex- 
pression for the threetangle has been observed in [4^ . 
giving an expression in terms of projectors) 



k3(x)l^ = ^^Ml,l^l^^M2,i^2^M3,l'3^-^M4,^'4®«A^I.(xkMl^M2f^K|X>(x|0-A'l^M2'^A|x) 
= 25KlAipil.i0K2A2M2l>2®K3A3M3l'3(xkKlO-K2O-K3lx)(xkAia'A2O-A3lx) 



(11) 

(12) 



Although this clearly favors the antilinear framework 
from a conceptual point of view, Eqs. (fTT|) . (fT2|) arc im- 
portant for practical purposes: they are a prescription 
for how the 3-tangle of a pure state could be directly 
measured in an experiment. 

An interesting conclusion from this approach is that 
it clarifies the relation between SL invariants and a pro- 
posal to use the projection operators P_ and P+ to as- 



sess entanglement in multipartite qubit systems |41| . To 
this end we note that the linear operator corresponding 
to CTyCfcf. Eq. (I6|)) equals precisely the operator P_ in 



Ref. l41[, which proves 5i(2,C) invariance of P_. On 
the other hand, P-^ is not SL invariant, which becomes 
also evident, when we transform it back to the anti-linear 
framework 



(V'|o(^|P+|V;)o|V^) ==7^,(V|ajV/)(V'*k,|^) ; 7^, = diag {1, 1, 0, 1} (13) 

I 

and take into account the results from Ref. [23|. As a sively from P_ in the way proposed in Ref. [4l| coincide 
consequence, an expression for qubits constructed exclu- with the well-known A^-tangle proposed by Wong and 



4 



Christensen (43| for an even number of qubits. Its use as 
a measure of multipartite entanglement has certain lim- 
itations p3l l43l . but an analytic convex roof extension 
is available [45|. This applies to every SL invariant con- 
structed in Ref. [4lj. A further limitation is that those 
quantities are polynomial local SU invariants with homo- 
geneous degree limited to 4. 

On the other hand, the mapping ([7|) of O = ay and re- 
sults from [s^l show how P_ can be used for constructing 
SL invariants of higher homogeneous polynomial degree. 
This procedure is entirely equivalent to employing full 
singlet contractions with the spinor metric tensor in the 
anti- linear framework (see e.g. [l3|); the latter method 
has been demonstrated to be considerably less efficient 
than using the pair of combs proposed in [23| . 

Experimental detection. ~ 

Now let us discuss the possibility of measuring a quan- 
tity such as the squared 3-tanglc in Eq. pT|) . As to ex- 
periments, one has to deal with mixed states p = (1 — 
£)|^')(^'| -f ep with some admixture p, where =0 
without loss of generality. This has two consequences: 
i) the 3-tangle typically will be reduced, Ti,{p) < T3(^'), 
and a) the quantity on the r.h.s. of Eq. ((TT|) looses its sig- 
nificance for growing e, as the experimentally measured 
correlation functions are expectation values for p rather 
than for (...) = tr(p. . .). In case of non-systematic 
errors, it should be reasonable, as a first approach, to 
assume this erroneous admixture to be completely un- 
correlatcd/unpolarized, i.e. proportional to the identity 
on the orthogonal complement of j'l'). For the discussion 
of unpolarized noise, we rewrite 



p=(l-e)|*)(vl'| + ^]l 



(14) 



where d is the dimension of the Hilbert space: for q qubits 
we have d — 2''. The identity matrix will only give a con- 
tribution for the identity operator, since any appearing 
Pauli matrix leads to trace zero. Then, the error of type 
a) is easily estimated as 



for the 2-tangle, i.e. the concurrence squared, and 



(15) 



rr = M^) - I (9r3(V) - 1 + m (16) 
for the 3-tangle, where we have defined 

M := 3(Al-)2 - 3{M,.)l + (M,,,)^ , 

(M,) := i^(^|a«|V>, 

3 ' 

(M_) := (V-lai^afai^)!^). 

This error applies when the tangle of the dominant pure 
state is computed in the laboratory using the expression 



of the pure state tangle in terms of spin expectation val- 
ues. 

Concerning the error type i), the admixture will typ- 
ically affect the multipartite entanglement in the state 
beyond the mere loss of weight e; this is generic when 
different entanglement classes coexist (see Ref. [i^) and 
will be referred to as entanglement interference. The lat- 
ter influence is reflected in the convex- roof extension, and 
this makes this second error considerably more complex. 
In order to estimate it, we discuss optimal decomposi- 
tions of weakly mixed states. 

Each decomposition vector of the density matrix ([M]) 
can be written as follows 



(17) 



where p = X)fcPfcl'/'fc)('/'fc| and J7 is a block of a unitary 
matrix, i.e. its column vectors Ufc := {Uo,k, Ui,k, . ■ .) are 
orthonormal; in particular = 1- For e = 

the optimal decomposition is trivially given by a sin- 
gle column vector with the only non-zero element being 
Uifi = 1. Continuity arguments for the decomposition 
vectors show that for £ ^ 1 a (5 =: Ae <g; 1 exists such 
that 



U 



( Pi 



\ 



\ 




(18) 



where X^jLi™ ~ 1 and without loss of generality 
/3 G H^. A will typically be of order 1. Then, 



m N — m 



As = 1 



i=l 1=1 

It is worth to notice that the matrix elements Si,j and 
7j J- will typically be of order one such that U itself will 
typically not be continuous. The convex roof of is its 
average in the optimal decomposition 

ra E(*'l'J'^)"'l('J'^l^'^22|$n(1'dS^''l'i>ni (19) 



There is one particular case to be distinguished: namely 
if m = 1 and Sij = 0. This case corresponds to a sit- 
uation, where the convex roof close to the state IvP) is 
affine on a range well beyond e. An example for this is 
a mixture of a GHZ state and an orthogonal W state 
as studied in [H, HJ, where A = ^jz^ ^ 1-17 with 
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pi = i + 3Yq-\/465 ~ 0.71. Then, the situation simpUfics 
a lot and a lower bound for the entanglement is simply 
given by r^lp] > r3(?/')(l — 2Ae). A drawback of this es- 
timate is however that in a generic experimental setup, 
it will be difficult to understand whether this scenario 
applies and what value A would have. 

For the estimation of T^[p] in a generic case, we define 
Xj := X]fc=i |f^j,fcP ^^'^ use the relation J^j^j = r — 
1. This situation corresponds to a situation where the 
convex roof close to j'l') is strictly convex. We need an 
estimate for the modulus of the states 



($,|cl>,) = |/3,.p(l-e) + 



k\ Pk 



(20) 



which can be given as 



< max{/3|(l - e),eXjPrnax} , 



(21) 



where pmin and Pmax are the minimal and maximal 
weights in j5; we will need the right inequality in what 
follows. In addition, we employ the very convenient 



estimate 



< 



3^y2edpma-xL , where wc have used that the largest eigcn- 

S]'^22 ^as 



-'(122 



value Amax of the tangle-operator ^ 
modulus 3. 

There appear two essentially different cases from (PT|) . 



If Aoi < 1 the 



discriminated by the ratio /3].ei := 
entanglement is dominated by p and the lower bound 
for its entanglement will be zero. When summing the 
condition /3j (1 — e) ^ sxjPmax over j we obtain (1 — 
s) ^ £Pmax(?' — 1) defining a rough threshold value for 
e beyond which the entanglement must be expected to 
have dropped to zero. For the analytically solvable case 
of GHZ-W-mixtures (r = 2, Pmax = 1) this threshold is 
£r = i . This rough estimate can only indicate an order of 
magnitude for such a threshold value. Indeed the analytic 
value for the threetangle to vanish is at e « 0.38 [i^. For 
unpolarizcd error of arbitrary rank {pmax = !/('' ~ l))i 
the threshold stays e = ^, and the correct value for r3 
to vanish is at e ~ 0.25 (see [i^ for specific r = 3 
states). 

Assuming now /^rci 3> 1 we obtain after some straight 
forward algebra that the lowest order deviation is 



£(t3[|*)] + 24(r - l)p 

max I • 



(22) 



For both the GHZ-W mixture and the fully unpolarizcd 
noise this deviation estimate would amount to e = 0.04. 
In all the above formulae, r denotes the rank of p. The 
details of the calculation will be presented elsewhere. The 
error clearly has a systematic contribution, namely due 
to the mere lack of weight of the original state, and a con- 
tribution due to entanglement interference as reflected in 
pccularities of the optimal decomposition. It must be 
noted that the interference term constitutes by far the 



most relevant term even for small e: assuming uncor- 
related noise in an experiment for three qubits we have 
r = 8 and the estimated maximal distructive entangle- 
ment interference accumulates to 24e. This would mean 
that en experimental error of only five percent could pos- 
sibly destroy the entanglement of a maximally entangle- 
ment GHZ state. In the cross-over region /Jioi ^ 1 we use 
< 2/3|(l — e) and the deviation up to first order 
in e is obtained as 



(23) 



For both the GHZ-W mixture and the fully unpolarizcd 
noise this deviation estimate would amount to e = 0.08. 
It is clearly seen that at least for weak mixtures to 
the GHZ state, both estimates are far too conservative. 
We must emphasize however, that a systematic analysis 
would be desirable in order to understand how conserva- 
tive the above estimates arc in general. 

In contrast to the error depicted in Eq. pB)) , the above 
formula attributes for the mixed state tangle of the state 
actually produced in the laboratory. It is important to 
note that there is no need to combine both errors to a to- 
tal error, since they correspond to different quantities to 
be measured. The first formula provides an experimen- 
tal check of whether the correct state has been produced 
(via its tangle characteristics) and a distillation proce- 
dure would be necessary to extract it from multiple copies 
in order to transmit it to quantum information process- 
ing. The second formulae rather address the question of 
how much entanglement there really is in the mixed state 
we have produced in the laboratory for direct submission 
to a quantum information process. 

Conclusions.- We present a one-to-one mapping which 
translates entanglement measures formulated within an 
antilinear formalism (23l . [23 | into expressions in terms of 
linear operators. These linear expressions have twice the 
multilinear degree of their more compact originals in the 
antilinear framework and thus constitutes an undue com- 
plication in a theoretical approach to entanglement quan- 
tification and classification which most naturally appears 
within the antilinear formulation. Thus, from a formal 
point of view our results favor an antilinear approach in 
order to obtain a deeper understanding of entanglement 
in multipartite qubit systems. This aspect is particularly 
important since for defining a proper entanglement order 
in multipartite systems each class of entanglement should 
be expected to have its own measure [11] . 

An important result of this work is an explicit expres- 
sion of the threetangle in terms of spin expectation val- 
ues. The transcription of multipartite entanglement mea- 
sures given as expectation values of antilinear operators 
in Ref.jl^, [IJ] is straightforward with the results pre- 
sented here, but it is not the scope of this work to list 
them explicitly. The invertible mapping ([7]) from the an- 
tilinear to the linear setting is not limited to applications 
for qubits. As soon as antilinear expressions for polyno- 
mial invariants of qudit systems are available, also these 
can be transcribed using the same technique. The linear 
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